Structural, electronic and magnetic properties of SrRuOs under epitaxial strain 
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Using density functional theory within the local spin density approximation, structural, electronic 
and magnetic properties of SRO are investigated. We examine the magnitude of the orthorhombic 
distortion in the ground state and also the effects of applying epitaxial constraints, whereby the 
influence of large (in the range of ±4%) in-plane strain resulting from coherent epitaxy, for both [001] 
and [110] oriented films, have been isolated and investigated. The overall pattern of the structural 
relaxations reveal coherent distortions of the oxygen octahedra network, which determine stability 
of the magnetic moment on the Ru ion. The structural and magnetic parameters exhibit substantial 
changes allowing us to discuss the role of symmetry and possibilities of magneto-structural tuning 
of SrRuOa-based thin film structures. 



I. INTRODUCTION 

Epitaxial thin film technology has the potential to 
change the functionality of conventional crystals in a 
way to achieve better performance and novel properties. 
Grown layer by layer in thin film geometries, the crys- 
tals still retain their bulk-like properties in two in-plane 
lateral dimensions, while being strongly modified along 
the direction of growth by finite-size effects, presence of 
interfaces and surfaces, and the misfit strain from epi- 
taxial lattice-matching. Additional interactions which 
emerge between consequent layers allow control, tuning 
and enhancing of their functional properties. Complex 
artificially-designed heterostructures with specific char- 
acteristics can be synthesized with atomic-level precision 
[!]• 

In many perovskite-oxide based heterostructures and 
superlattices, the metallic ferromagnetic perovskite 
SrRuOa plays a particularly key role, often serving as an 
electrode material that allows for better integration and 
facilitates electrical measurements [2, 3]. However, recent 
studies of epitaxial thin films [4] suggest that SrRuOs 
may have novel magnetostructural properties in ultra- 
thin film form. By tuning film thickness and in-plane 
strain (through choice of substrate) very different prop- 
erties may emerge as compared with the bulk. Simi- 
lar approaches have been exploited in recent experimen- 
tal and theoretical investigations showing modifications 
of the ferroelectric epitaxial thin films with strongly en- 
hanced performance [5-9]. However, to our knowledge, 
the impact of strain on ferromagnetic perovskite films has 
not yet been studied. 

Bulk SrRuOs undergoes a series of phase transforma- 
tions with decreasing temperature (for more details see 
Ref. [10, 11]). It is a metallic ferromagnet with a Curie 
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temperature about 160 K and a magnetic moment of 1.1 
/is per formula unit [12]. The ideal cubic perovskite 
structure with cubic symmetry PmSm is stable above 
950 K. From this temperature and down to 820 K it ap- 
pears in the tetragonal I A/ mem structure. Below 850 K 
and at standard pressures, SrRuOa has an orthorhombi- 
cally distorted perovskite structure (space group Pbnm 
[41] with 20 atoms per cell [13]). The structure is the 
GdFeOa / CaTiOa -type, which can be understood by first 
doubling the 5-atom cell along [110] with both oxygen oc- 
tahedra being rotated in the same way around [001] (z- 
axis), and then doubling this 10-atom system along [001] 
with two oxygen octahedra connected along the ^-axis, 
being further tilted in opposite directions. This structure 
is shown in Fig 1. In the Glazer notations, the tilting of 
SrRuOa is described by {a~a~c'^) [14, 15]. One of our 
main points below will be that in SrRuOa, the rotation 
and tilting is accompanied by a substantial deformation 
of the oxygen octahedra. 

Previous first-principles studies of SrRuOs examined 
the electronic and magnetic properties of the bulk low- 
temperature Pbnm structure [16-18]. Their results 
largely agree with experimental observations and explain 
basic features of the ground state in the absence of strain. 
In this paper, we use first-principles calculations to iso- 
late and investigate the effect of epitaxial strain on the 
structure of SrRuOs. We explore SrRuOa with this con- 
straint through examination of the orthorhombic distor- 
tion, magnetic and electronic properties of the Pbnm 
structure, elucidating the diiferences between bulk and 
epitaxially strained thin films in the [GDI] and [110] ori- 
entations. 

In section II we give a short description of our methods 
used in this work. Section III presents our most impor- 
tant results, which we split into two parts. In part (A) 
we concentrate on the structural parameters of a strain- 
free structure and highlight a confusion in the literature 
about the lattice parameters of SrRuOs. Part (B) deals 
with [001] and [110] oriented structures of SrRuOs grown 
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FIG. 1: (Top picture) The structure of the oxygen network in 
a Pm3m perovskite structure. The oxygen octahedra are per- 
fectly regular. (Bottom picture) The structure of orthorhom- 
bic Phnm perovskite structure with tilting (rotation of the 
octahedral cages). The oxygen octahedra in this case are not 
necessary regular. We consider this aspect in section III of 
this paper. 



with tensile and compressive misfit strain. Section IV 
presents discussions of three selected topics which follow 
from our results and deserve more attention. Conclusions 
are given in Section V. 



II. METHODS 

All calculations are performed with the Vienna ah- 
initio Simulations Package(VASP)[19, 20], with the 
Ccpcrley- Alder parametrization of the local spin-density 
approximation (LSDA) and projector-augmented wave 



potentials [21, 22]. The calculations arc performed us- 
ing a plane-wave energy cutoff of 500 eV and the conver- 
gence was checked with a cutoff of 600 eV. For calcula- 
tions with the Phnm unit cell of SrRuOa, which contains 
4 formula units (f.u.) or 20 atoms per cell, a 12 x 12 x 10 
Monkhorst-Pack fc-point mesh is found to yield conver- 
gence of all properties computed here and is used in all 
self-consistent calculations. For a simple cubic perovskite 
cell with 5 atoms we used a 12 x 12 x 12 fc-point sampling. 

In order to isolate the effects of epitaxial strain, we per- 
form bulk calculations of the Phnm ground-state struc- 
ture with lattice vectors constrained to match a hypo- 
thetical substrate. Specifically, for a series of values of 
misfit strain, we fix the in-plane {xy plane) lattice vec- 
tors, allowing the c lattice parameter (perpendicular to 
the substrate) and the internal ionic coordinates to fully 
relax within the symmetry of the Phnm space group. 
The ions were free to move until the Hellmann-Feynman 
forces were less than 4 meV/A. 

Two different angles characterize the degree of rota- 
tion and tilting of the oxygen octahedra in SrRuOs. The 
rotation and tilting represent basically the same type of 
distortion, but with respect to different directions. It is 
conventional to define the rotation along the z axis sep- 
arately from those along the x and y axes, referred to as 
tiltings. For what follows, we define the dihedral angle 
<d to be along the RU1-O2-RU3 bond shown in Fig. 2, and 
the tilting angle is given by (180-9)/2. Similarly, the 




FIG. 2: The Phnm crystal structure of SrRuOs (GdFeOa 
/CaTiOa -type). The unit cell consists of four formula units 
of the ideal cubic perovskite structure. The Ru atoms possess 
relatively large magnetic moments and occupy high-symmetry 
positions in the orthorhombic unit cell. The atoms of oxygen 
and Sr are displaced from their high-symmetry positions due 
to the tilting. The numbers attached to certain atoms help 
us to define the angle of tilting as discussed in the section 
"Methods" . 
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FIG. 3: (Top picture) The top view of the structure of SrRuOa 
showing the distortion which we refer to as a rotation. The 
corner-connected oxygen octahedra rotate in opposite direc- 
tions around the [001] direction. This distortion changes the 
a and b parameters of the structure equally. (Bottom pic- 
ture) The tilting of the SrRuOs structure shown projected 
along the [110] direction (or 6 axis of the Phnm structure). 
In this view one can see the folding of the [001] plane that is 
responsible for the shrinking of the Pbnm structure along the 
a axis. 



rotation angle is defined through the relation (90-r2)/2, 
where VI is an angle shown in Figures 2 and 3. The rota- 
tion and tilting commute with each other and uniquely 
describe the network of oxygen octahedra, assuming the 
octahedrons are regular. We will return to the validity 
of this assumption later in section III [42]. 



III. RESULTS 
A. Strain-free structure of SrRuOs 

Starting from the experimental space group Phnm, we 
first obtain the ground state structure of SrRuOa by re- 
laxing all seven free internal structural parameters (in- 
ternal atomic coordinates) and three lattice parameters. 
Our results, along with a comparison with the experi- 
mental values for different perovskite systems, are shown 
in Tables I and II. The calculated parameters are in good 
agreement with experiment. The slight underestimate of 
the lattice constant is typical for LSDA. We calculate 
the FM state in Phnm structure to have an energy 186 
mcV/f.u. lower than that of the cubic perovskite FM sys- 
tem, which was optimized within the Pm3m symmetry 
(see Table III). A previous first-principles study (Rcf. 16) 
computed a comparable energy difference, 140 meV/f.u., 
using structures with experimental lattice parameters. 

Experimental values of the magnetic moment in 
SrRuOa range between 1.1 and 1.6 /xs/f.u.; the spread 
in values is expected given the difficulty of making 
single-domain samples and their large magnetocrystalline 
anisotropy [12, 23]. According to Ref. 18, only approxi- 
mately 60% of the total magnetic moment is Ru derived, 
while the remaining 40% is distributed among oxygen 
sublattice sites. A similar spread of the values of the 
magnetic moment in SrRuOa has been reported from 
theoretical investigations. Usually they are obtained us- 
ing the LSDA and experimental lattice parameters for 
the orthorhombic Pbnm structure. Such calculations 
reported magnetic moments ranging between 0.97 and 
1.96 /is/f-u. [16, 17, 24, 25] 

With an intermediate degree of relaxation, keeping the 
experimental lattice parameters fixed, but allowing the 
Wyckoff positions to change, we obtain 1.2 /is/f.u. Af- 
ter a complete relaxation of the structure, our magnetic 
moment calculated with LSDA is 0.98 /is/f.u., smaller 
than that obtained with experimental lattice parameters, 
as would be expected with the smaller LSDA lattice pa- 
rameters. A summary of all calculated lattice parameters 
appears in Table II. The LSDA slightly overestimates 
bond strength and therefore favors stronger hybridiza- 
tion of the Ru-d and 0-p states, giving smaller lattice 
parameters and consequently smaller magnetic moments. 

In addition to the ferromagnetic (FM) configuration, 
we optimize the bulk structural parameters of SRO with 
different initial magnetic orderings, including the non- 
magnetic (NM) phase and three different antiferromag- 
netic (AFM) spin arrangements. Several local minima 
in the total energy are found. Among these minima, 
the FM configuration is computed to be the most stable 
state, in agreement with experiment. The A- type AFM 
and C-type AFM, and NM configurations are found to 
be 6.77meV/f.u., 6.44meV/f.u., and 8.15meV/f.u. higher 
in energy than the FM ground state, respectively. The 
G-type AFM configuration is found to be unstable: when 
relaxed, it reverts to the NM minimum. In the remaining 
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TABLE I: Experimental structural parameters for several perovskite systems with Pbnm symmetry taken from Ref. 13, 26-30. 
We focus in this list on the relations of the a, b and c lattice parameters. In most of the systems a < 6 is satisfied due to the 
tilting of the oxygen octahedra, whereby only the a vector shrinks, while the b vector remains unchanged. From this table, one 
can see that for SrRuOs the ratio a < 6 is not satisfied. 





a(A) 


6(A) 


c(A) 


Wyckoff position 


X 


y 


z 


SrRuOs 


5.5670 


5.5304 


7.8446 


Sr (4c) 


-0.0027 


0.0157 


0.25 










Ru (4b) 


0.5 


0.0 


0.0 










O (8d) 


0.7248 


0.2764 


0.0278 










O (4c) 


0.0532 


0.4966 


0.25 


GdFeOa 


5.3510 


5.6125 


7.6711 


Gd (4c) 


-0.0015 


0.0626 


0.25 










Fe (4b) 


0.5 


0.0 


0.0 










(8d) 


0.6966 


0.3011 


0.0518 










(4c) 


0.1009 


0.4669 


0.25 


GdAlOs 


5.2537 


5.3030 


7.4434 


Gd (4c) 


-0.0079 


0.0376 


0.25 










Al (4b) 


0.5 


0.0 


0.0 










(8d) 


0.7147 


0.2855 


0.0387 










(4c) 


0.0724 


0.4863 


0.25 


CaTiOa 


5.3804 


5.4422 


7.6417 


Ca (4c) 


-0.0065 


0.0349 


0.25 










Ti (4b) 


0.5 


0.0 


0.0 










(8d) 


0.7111 


0.2884 


0.0372 










(4c) 


0.0707 


0.4842 


0.25 


CaRuOa 


5.3408 


5.5311 


7.6460 


Ca (4c) 


-0.0150 


0.0560 


0.25 










Ru (4b) 


0.5 


0.0 


0.0 










(8d) 


0.6930 


0.2970 


0.0530 










(4c) 


0.0910 


0.4670 


0.25 


GdScOs 


5.4862 


5.7499 


7.9345 


Gd (4c) 


-0.0163 


0.0594 


0.25 










Sc (4b) 


0.5 


0.0 


0.0 










(8d) 


0.6931 


0.3007 


0.0556 










(4c) 


0.1183 


0.4465 


0.25 


DyScOa 


5.4400 


5.7130 


7.887 


Dy (4c) 


-0.0172 


0.0607 


0.25 










Sc (4b) 


0.5 


0.0 


0.0 










(8d) 


0.6926 


0.3040 


0.0608 










(4c) 


0.1196 


0.445 


0.25 



TABLE II: Calculated structural parameters for several perovskite systems with Pbnm symmetry. Theoretical results also show 
the distortion of SrRuOs for which the orthorhombic distortion is opposite to other perovskites. We introduce parameters, 
Aa and Ab, which characterize shrinking of the a and b lattice parameters, respectively. The parameter a shrinks due to the 
tilting, while b shrinks due to the rectangular distortion of the horizontal middle plane of the oxygen octahedra. 





a(A) 


KA) 


c(A) 


Wyckoff pos. 


X 


y 


z 


Tilting 


Aa 


Ab 


SrRuOs 


5.5031 


5.4828 


7.7546 


Sr (4c) 


-0.0050 


0.0296 


0.25 


10.21° 


1.58 


1.75 










Ru (4b) 


0.5 


0.0 


0.0 
















(8d) 


0.7165 


0.2834 


0.0336 
















(4c) 


0.0647 


0.4941 


0.25 








CaRuOs 


5.2090 


5.5297 


7.5512 


Ca (4c) 


-0.0211 


0.0645 


0.25 


16.79° 


4.26 


2.15 










Ru (4b) 


0.5 


0.0 


0.0 
















O (8d) 


0.6929 


0.2996 


0.0519 
















O (4c) 


0.1034 


0.4665 


0.25 








CaTiOs 


5.2900 


5.4007 


7.5334 


Ca (4c) 


-0.0099 


0.0468 


0.25 


13.26° 


2.67 


0.75 










Ti (4b) 


0.5 


0.0 


0.0 
















O (8d) 


0.7065 


0.2926 


0.0425 
















O (4c) 


0.0811 


0.4790 


0.25 
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FIG. 4: This picture illustrates schematically the fact that 
as a result of tilting the Pbnm structure is necessarily or- 
thorhombic. This appears from simple geometrical consider- 
ations when tiltings a~a~ have the effect of folding of the 
Ru-0 (001) atomic planes. If the folding develops along the 
a axis, while another lattice vector 6 remains unchanged, we 
have a < b. The importance of this point can be seen when the 
Pbnm structure is subject to the epitaxial constraints enforc- 
ing a — b. The rotation around the z axis is different because 
it changes the a and b parameters of the Pbnm structure in 
equal proportions. So, the a = b constraint requires the mid- 
dle plane of the oxygen octahedra to shrink along the b axis. 
In order to characterize this epitaxially constrained structure 
we introduce the notation e-Pbnm. The angles shown de- 
fine degrees of tilting and rotation as discussed in the section 
"Methods" . 



of this paper, we focus on the FM state and its compar- 
ison with the NM state. 

At this point, we turn to a discussion of the lattice pa- 
rameters of the Pbnm structure. From Tables I and II, we 
can see that most of the Pbnm GdFeOs-type structures 
have a < b < c. Interestingly, however, in the case of 
SrRuOa this relation is not satisfied. SrRuOs is the only 
compound we know of for which this is the case. This 
has led to a confusion in the literature where one can find 
reports both of a > 6 (see Refs.ll, 13, 31-33) and a < b 
(see Refs.23, 24, 34, 35); in the latter case parameters are 
apparently switched to be consistent with overall trend in 
the GdFeOs-type structures. To clarify the unique situa- 
tion in SrRuOa other geometrical factors making up the 
orthorhombic distortions of the Pbnm structure need to 
be considered, as we now explain. 

Figure 3, together with a sketch in Fig. 4, illustrates the 
relationship between octahedral tilting and rotation and 
the orientations of the interatomic bonds. The in-plane 
lattice parameters of the Pbnm structure are not equal, 
that is a 7^ &, as a consequence of the tilting; the rotation 
on the other hand does not change the ratio of a to b. 
For arrangements with Pbnm symmetry, if the octahe- 
dra were regular, tilting would require a < b, and not 
the other way around. This simple geometrical picture 
explains the orthorhombic shape of most of the known 
perovskite systems listed in Table I, with the notable ex- 
ception of SrRuOs. 

Our theoretical calculations confirm that for SrRuOa, 
a > b (see Tab. II), while for the systems like CaRuOs 



and CaTiOs, a < b. An analysis of the results of our cal- 
culations shows that the reason for the apparent anomaly 
in SrRuOs is that the oxygen octahedra are not regular: 
the "horizontal" middle plane is not square, but rather 
is rectangular. The rectangular distortion occurs along 
the b axis partially compensating the orthorhombic dis- 
tortion originating with the tilting. Thus, two different 
contributions determine the orthorhombic distortion (the 
a/b ratio) in the Pbnm structure. 

In Table II we summarize the impact of distortions 
coming from each contribution. The anomalous or- 
thorhombic distortion in SrRuOs is derived from geo- 
metrical considerations and those parameters which we 
already have: a, b, c, tilting and rotation angles. If a°'^^ 
and 6°'^' describe the equatorial plane of the octahedra, 
then the lattice parameters of Pbnm can be expressed as 
follows: 



a = 2a°'^' cos (rotation) cos (tilting) ; 
b = 26°'^* cos (rotation). 



(1) 



Consequently, the orthorhombic distortion can be ex- 
pressed as: 



5oct 



cos (tilting). 



(2) 



These simple relations show that the internal distor- 
tions can be readily determined from fundamental pa- 
rameters of the Pbnm structure. However, it is impor- 
tant to note that the oxygen octahedra are not necessarily 
regular, i.e. a°'^^ and 6°*^' are free parameters unique for 
each material. 



B. Epitaxially constrained structure: e-Pbnm 

We now consider the structure and properties of 
epitaxially strained Pbnm SrRuOa. Experimentally, 
SrRuOs films can grow with either [001] or [110] ori- 
entations on various substrates like LaAlOs, SrTiOa, 
DyScOa or GdScOs [33, 34] (see Ref. 35 for good illus- 
trations of the two geometries). The [001] case is a rel- 
atively straightforward extension of previously studied 
perovskite multi-layers [36, 37], except that additional 
degrees of freedom must be considered to account for the 
tilting and rotation. In the [110] orientation, the situa- 
tion is quite different, as we discuss in detail below. 



1. Geometries of the [001] and [110] oriented films 

For the [001] oriented films, the a and b parameters 
of the orthorhombic Pbnm lattice are constrained to be 
equal (in order to match the square lattice of the sub- 
strate). We refer to these constrained structures, with 
all other structural parameters relaxed, as e-Pbnm[Q01\, 
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[110] 




FIG. 5: Schematic presentation of tlic epitaxial structure of 
SrRuOa grown in the [110] orientation. The grey parallelo- 
grams represent the equatorial planes of the oxygen octahedra 
in the (001) plane. This geometry implies that at zero in-plane 
strain the structure still have the orthorhombic symmetry, but 
its [110] direction is not normal to the substrate. Due to the 
orthorhombic symmetry, there is an angle 7 defined by eq.(3). 
Applying in-plane strain changes this angle, but at the same 
time the orthorhombic symmetry is broken. SrRuOa grown 
with the [110] orientation under strain is monoclinic. 
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FIG. 6: The comparison of the elastic energies of FM e- 
P6nm[001] SrRuOs and NM e-P6nm[001] CaTiOa as a func- 
tion of the applied misfit strain. The in-plane lattice param- 
eters were fixed by the value of the misfit strain, while all 
other parameters of the structures were relaxed to the zero- 
force states. In addition we show here the total energy of the 
simple PA/mmm perovskite SrRuOs when rotation and tilt- 
ing are not allowed. This allows us to conclude that the tilting 
does not change the elasticity of the structure with respect to 
the misfit strain. 



with the [001] orientation. The diagonal is fixed 

and the lattice parameter c is the in-plane parameter. 
The out-of-plane lattice vector is not perpendicular to 
the surface (see Figure 5), because its angle with the 
normal to the surface is given by the difference between 
parameters of Pbnm, a and b: 



where "e" indicates "epitaxial" . Misfit strain is mea- 
sured relative to the lowest energy e-P6nm[001] struc- 
ture, the computed parameters being given in Table III. 
Note that the zero misfit strain structure is not the same 
as the strain-free structure discussed in the previous sub- 
section, as the latter is not compatible with the epitaxial 
constraints. For the e-Pbnm structure we have Aa = Ab 
(from Table II), i.e. two kinds of distortions compensate 
each other, as they are forced to by the imposed condi- 
tions. 

In the [110] orientation, we consider a general case 
when the SrRuOs film is grown on an orthorhombic sub- 
strate. We assume that our substrate has in-plane non- 
equal lattice parameters, but they are kept orthogonal. 
If we made the in-plane parameters equal, this would 
correspond, for example, to a cubic lattice of SrTiOs 
with a compressive lattice mismatch of about 0.59%. But 
making the in-plane lattice rectangular brings us close to 
the situation with GdScOa or DyScOa substrates, whose 
structures match nicely with [110] oriented structure of 
SrRuOa (we refer to it as e-P6rim[110], although we show 
that it may have lower symmetry). In this case SrRuOs 
is subject to substantially different constraints compared 



7 = 2 arctan J =2 arctan ^ cos (tilting) J (3) 

Although the diagonal ^/a^ + b'^ is fixed on the sub- 
strate, the ratio a/b is free to change. From eq.l we 
see that b/a = (&°'^*/a°'^*) cos (tilting). At zero in-plane 
strain the [110]-oriented SrRuOa retains the bulk struc- 
ture if it is grown on SrTiOs, but not in the case of [001] 
oriented GdScOa or DyScOa substrates. 



2. Structural changes as a function of strain 

As the epitaxial constraint is varied, the shape of the 
cell is changed, which increases its internal energy. In 
Fig. 6, we plot the elastic energy of ferromagnetic SrRuOa 
for misfit strains ranging from -4.5% to +3.5%. Sur- 
prisingly, the energetic penalties for the in-plane distor- 
tion are nearly the same for both e-Pbnm and PA/mmm 
structures, though the latter has no internal degrees of 
freedom to relax and reduce the energy. Moreover, com- 
paring these curves with another perovskite, CaTiOs, we 
see that their behavior is strikingly similar. Changes of 
the volume are shown in Fig. 7. 
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TABLE III: Computed structural parameters of SrRuOa in the cubic perovskite (space group Pm3m) and the epitaxially 
constrained structure e-Pbnm[001] in the zero strain state. Results for both FM and NM configurations are listed in order to 
show that they do not have significant differences (later on, it will be shown that applying epitaxial strain yields important 
differences between the FM and NM cases). The Wyckoff positions for all structures refer to space group Pbnm. 





a(A) 


fe(A) 


c(A) 


Wyckoff pos. 


X 


y 


z 


Tihing 


Aa 


Afe 


e-Pbnm{FM) 


5.4929 


5.4929 


7.7514 


Sr (4c) 


-0.0055 


0.0303 


0.25 


10.48° 


1.67% 


1.67% 










Ru (4b) 


0.5 


0.0 


0.0 
















(8d) 


0.7160 


0.2837 


0.0336 
















O (4c) 


0.0650 


0.4937 


0.25 








e-P6nm(NM) 


5.4929 


5.4929 


7.7440 


Sr (4c) 


-0.0051 


0.0304 


0.25 


10.48° 


1.67% 


1.69% 










Ru (4b) 


0.5 


0.0 


0.0 
















O (8d) 


0.7155 


0.2842 


0.0345 
















(4c) 


0.0650 


0.4945 


0.25 








Cubic(FM) 


5.5070 


5.5070 


7.7880 


Sr (4c) 


0.0 


0.0 


0.25 


0.0 


0.0 


0.0 










Ru (4b) 


0.5 


0.0 


0.0 
















(8d) 


0.75 


0.25 


0.0 
















O (4c) 


0.0 


0.5 


0.25 









The degree of tilting and rotation of the oxygen octahe- 
dra reflects the impact of epitaxial strain. Figure 9 shows 
the evolution of these angles for both [001] and [110] ori- 
entations. For [110] the two angles decrease together. 
On the other hand, for the [001] orientation, the angles 
behave oppositely: for increasing compressive or tensile 
strain, the rotation or tilting become more pronounced, 
respectively. However, we find that the angles exhibit 
surprising behavior under tensile strain in the NM state 
(discussed further below). 



IV. DISCUSSIONS 



From the results reported above, we selected three par- 
ticular points for further discussion. Our focus here is on 
the magnetostructural coupling in SrRuOs and its de- 
pendence on the film orientation. 



In addition to the orientation of the oxygen octahcdra, 
we observe their shape is also changing. This is especially 
important because we have to understand the behavior of 
the angles in the NM case shown in Figure 9. The split- 
ting of two Ru-Oz curves in Figure 10 can be explained by 
the fact that the NM configuration allows for additional 
tetragonal contraction of the oxygen octahedra along the 
z axis. This contraction manifests itself in the trend in 
angles shown in Fig. 9. Connection of this effect with the 
magnetic properties of SrRuOa will be discussed later in 
Section IV. 



Applying strain induces a change in volume. In the 
[110] orientation, the volume can be modified by rotating 
the angle 7 or by expanding (or contracting) the out- 
of-plane lattice parameter. The latter can be achieved 
only if the symmetry of the oxygen octahedra is broken. 
Namely, the middle plane Ru-0 bonds would be required 
to have different length; in the case of [001] oriented films 
they are the same (see Fig. 10). This additional lowering 
of the symmetry alters the electronic environment of the 
Ru ion leading to a drop of its magnetic moment, as 
shown in Fig. 8. The [110] oriented structure of SrRuOa 
becomes monoclinic [43] (see Fig. 5). 



A. Magnetoelastic coupling 

Fig. 7 exhibits some nonlincarity in the volume of 
the PA/mmm with misfit strain, in contrast with the 
e-P6nm[001] structure. To understand the behavior of 
PA/mmm it is useful to compare Figures 7 and 8. Both 
curves show inflection points that in Fig. 7 coincide with 
abrupt drops in magnetic moment. Evidently, the mag- 
netic moment, via the spin configuration of the Ru ion, 
stabilizes a symmetrical non-deformed shape of the oxy- 
gen octahedra. Concurrently, the in-planc strain leads 
to the out-of-plane lattice expansion (contraction) under 
the compressive (tensile) strain. At some critical strain 
the contraction becomes more favorable than the mag- 
netic orbital ordering and the lattice distorts. On the 
other hand, in the e-Pbnm structure, the tilting and ro- 
tation preserve the shape of octahedra and the magnetic 
orbital ordering survives. 

The behavior described here suggests a strong magne- 
tostructural interaction, which may potentially explain 
past experiments showing the anomalous thermal expan- 
sion of SrRuOa [28]. We reserve a more complete inves- 
tigation for future work. 
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FIG. 7: Lattice parameter c and the volume of the [001] and 
[110] oriented perovskite structures of SrRuOs under epitaxial 
strain. (Top) The epitaxially constrained e-Pbnm structure. 
Both volume change and tetragonal distortion cost energy, 
whereby the structure changes these degrees of freedom si- 
multaneously. Both [001] and [110] oriented films show simi- 
lar changes in the volume as a function of the strain. 
(Bottom) Non-tilted [001] oriented simple perovskite 
PA/mmm. Both curves show remarkable difference of their 
behavior as compared to the e-Pbnm structure. A compari- 
son to the Fig. 8 shows that the distortions of the P4:/mmm 
are strongly influenced by the magnetoelastic coupling. 



B. Anomaly of tilting and FM-NM transition 

In the previous section we saw the effect of coupling 
between the distortions of the oxygen octahedra and the 
stabihty of the magnetic moment on Ru in the PA/mmm 
symmetry. This allows us to understand what we see in 
the case of c-P6rim[001] structure. In Figure 10, we see 
an extra contraction of the Ru-Oz bond, which occurs 
in the NM phase under tensile strain. In contrast, in 
the FM phase the oxygen octahedra keep their shape by 
monotonically increasing the tilting angle. 

This difference in the behavior of the magnetic and 
non-magnetic phases can be explained from a simple spin 
configuration argument known for SrRuOa [38]. If the 
structure is not distorted, four 4c? electrons of Ru'*'*' oc- 
cupy 3- fold degenerate t2g orbitals leading to (T^icc i^2g)' 
whereby two spins are not compensated and the Ru ion 
has a large magnetic moment. If the oxygen octahedra 



0.5 



O-o e-Pbnm [001] 
P4/mmm [001] 
a-<2 e-Pbnm [110] 




-2 
In-plane strain (%) 

FIG. 8: The total magnetic moment of SrRuOs as a function 
of strain. In the case when tilting is not allowed (PA/mmm), 
the deformation of the Ru-O octahedra, occurring due to the 
strain, becomes so strong that the total magnetic moment 
drops down to zero. Redistribution of the electron orbitals 
changes the magnetic interactions. Beyond critical values of 
the strain, SrRuOs becomes non-magnetic. When the tilting 
is present {Pbnm structure) the magnetic moment is retained 
even at large in-plane strain. In the [110] oriented film the 
magnetic moment is significantly smaller, which we attribute 
to the geometry of the oxygen cages and the corresponding 
changes in the electronic structure (see DOS in Fig. 12). 



are contracted along some of the Ru-O bonds, the t2g or- 
bitals split producing a non-magnetic spin configuration 
forRu, "itlgd^g- 

Such a transition is possible and can have signifi- 
cant implications for the magnetostructural tuning of the 
SrRuOa based hcterostructures. One could change the 
magnetic ground state of SrRuOs by using uniaxial com- 
pression. Starting from the FM state, one can apply a 
uniaxial stress along the c axis. The NM state which 
has a shorter lattice parameter c can become the ground 
state. Some doping of the SrRuOs structure with larger 
atoms, like Ba, should make tilting weaker and therefore 
favor the NM state, i.e. reduce the amount of pressure 
required to transform the structure. We will discuss this 
issue in more detail with possible applications in our next 
paper. 



C. Orientation dependence 

We now comment on the differences between [001]- and 
[110]-oricnted films. For the [001]-oriented films, the oxy- 
gen octahedra are free to rotate and tilt adopting to the 
changing lattice parameters. In contrast, in the [110] 
orientation both angles are fixed to the parameters of 
substrate. Moreover, we stressed above that the [110] 
oriented SrRuOa becomes monoclinic. This difference 
has consequences for the calculated electronic structure. 
We show the electronic DOS projected on the Ru sites 
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FIG. 9: Tilting and rotation angles of the oxygen octahe- 
dra as a function of strain. (Top) The e-Pbnm[001] structure 
shows a sudden change in the behavior of the tilting and ro- 
tation angles in the NM system, which turns out to be due 
to the additional contraction of the oxygen octahedra along 
the z axis (see Fig. 10). Both in the FM and NM cases un- 
der tensile strain the rotation angles stay close to a value of 
10 degrees. (Bottom) In the e-P6nm [110] orientation, both 
angles are fixed by the substrate constraint. 



for four different cases. These four plots arc intended to 
illustrate how the electronic environment of Ru changes 
when distortions of the oxygen cage network occur. 

The additional lowering of the symmetry in the [110] 
orientation leads to a splitting of the t2g bands associated 
with Ru d level, and a smaller magnetic moment (Fig. 8). 
The gap opens at about 0.8 eV below for the spin-up 
states, and 0.4 eV below E^;^ for the spin-down component 
(Fig. 12(Bottom)). 

From examination of the electronic structure in Fig. 
11 and 12 four additional features emerge. We summa- 
rize them as follows: First, for Pm3m, there are no dis- 
tortions of the oxygen octahedra, and thus there is no 
gap in DOS (see Fig. ll(Top)); second, for Pbnm the 
octahedra are rotated, thus the t2g and Cg bands start 
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FIG. 10: Bond lengths of Ru-0 inside the oxygen octahedral 
cages as functions of in-plane strain. (Top) The [001] oriented 
film with FM and NM configurations. Under the tensile strain 
the NM and FM curves split, which is due to the additional 
contraction of the oxygen octahedra along the z axis in the 
NM case, which is reflected in smaller tilting angles shown in 
Fig. 9. (Bottom) In the [110] orientation, the substrate has an 
orthorhombic shape, therefore the in-plane parameters of the 
SrRuOs film are distorted. In order to be consistent we use 
the notations of the [001] orientation. Therefore, the x axis 
in this case is out of plane. 



to "feel" each other and repel, and as a consequence, 
there is a gap just above Ep in the DOS (see Figure 
11 (Bottom)); for e-Pfenm [001] under tensile strain, the 
rotation is reduced and the i2g{xy) and eg(x^ — y^) re- 
cover their symmetry, whereby there repulsion becomes 
weaker, reducing the gap above Ep (see Fig. 12(Top)). 
Finally, for e-P6nm[110], this orientation requires further 
lowering of the symmetry. The crystal field splits the t2g 
bands and opens an additional gap 0.8 eV (0.4 eV for 
spin-down) below the Ep. 
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FIG. 11; Ru site projected partial DOS of SrRuOa. The 
electronic states with Eg symmetry are filled with gray color. 
Different orientations of the t2g and e^ orbitals are plotted to- 
gether without distinction in order to focus on the gaps which 
open when distortions are applied to the structure. (Top) 
DOS of the simple perovskite structure with PmSm symme- 
try. There is no gap above the Fermi level. (Bottom) DOS 
of the strain free bulk Phnm structure. There is a gap above 
the Fermi level, which opens, as we discuss in the text, due 
to the rotations of the oxygen octahedra. In Fig. 12 it will 
be shown that in the [110] oriented structure, one more gap 
opens. 



V. CONCLUSIONS 

In this paper, we studied properties of SrRuOa under 
conditions of epitaxial strain and different thin film orien- 
tations and considered how the changes that take place 
in SrRuOa may be used to tune properties of SrRuOa 
and SrRuOa based heterostructures. The results can be 
summarized as follows: 

We explained the mechanism leading to the or- 
thorhombic distortion of Phnm, namely, the tilting and 
the contraction of the oxygen octahedra. The "anoma- 
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FIG. 12: Ru site projected partial DOS of SrRuOs in the 
vicinity of Ef- The electronic states with eg symmetry are 
filled with gray color. Different orientations of the t2g and 
eg orbitals are plotted together without distinction in order 
to focus on the gaps which open when distortions are applied 
to the structure. (Top) DOS for the e-P6nm[001] structure 
of SrRuOa under tensile strain of 2.5%. The tensile strain 
reduces in-plane distortion of the oxygen octahedra, whereby 
the gap becomes smaller. (Bottom) DOS for the e-Pfenm [110] 
structure. The strain is now also tensile, but 1.94%. There 
is another gap which opens about 0.8 eV (0.4 eV for spin- 
down) below the Fermi level. It happens because the [110] 
orientation required additional distortion of the oxygen cages, 
forcing the t2g electrons of Ru to split. 



lous" orthorhombic shape of SrRuOa can be easily un- 
derstood from simple geometrical considerations. 

Two different orientations of the SrRuOa epitaxial 
films, e-P6rim[001] and c-Phnm[\\Q\, exhibit significant 
differences in their structural properties. The latter, un- 
less there is no in-plane strain, has monoclinic symmetry, 
while the [001] oriented film remains orthorhombic at any 
reasonable value of the strain. 

SrRuOs exhibits significant magnetostructural cou- 
pling. Tensile strain reveals that in the non-magnetic 
state there is a distortion of Ru-Oz bonds which can be 
attributed to the change of Ru spin configuration, ex- 
plaining the fact that this distortion does not appear in 
the magnetic state. The tilting and rotation help to pre- 
serve the shape of the octahedra when epitaxial strain 
is applied, whereby the magnetic orbital ordering is pre- 
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served. By suppressing the tilting under strain one could 
obtain a non-magnetic state of SrRuOa. 
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